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solution for the deflections, moments, and reactions rectangular plate 
under transverse load, supported the corners, and with without supporting 
beams the panel boundaries, presented this paper. This solution can 
classified unified one, since adaptable most the specific problems 
for which particular solutions have been developed. Provision made for 
other classes particular solutions for additional cases. The solutions for 
specific problems, expressed general form, differ only the arbitrary 
constants appearing the general form. These constants are evaluated 
means operational formulas, with dependence parameters that are defined 
the ratios the flexural and torsional rigidities the beams the flexural 
rigidity the plate, and the proportions the plate. illustrate the 
procedure, the arbitrary constants are obtained for rectangular plate under 
uniform load that simply supported beams finite flexural rigidity 
four sides. particular value the parameter relative flexural rigidity for 
which the expressions for moment reduce simple exact formulas also exists. 
Moment and reaction coefficients are evaluated positions along the axes 
symmetry and the boundaries the simply supported plate for several values 
the parameter relative flexural rigidity. Experimental verification 
obtained use introductory tests panels supported beams 
arbitrarily reduced flexural rigidity. 


INTRODUCTION 


Notation.—The letter symbols introduced this paper are defined where 
they first appear, the text illustration, and are assembled alphabeti- 
cally Appendix for convenience reference. 


comments are invited for publication; the last discussson should submitted 
November 1953. 
Analyst, Loewy Construction Co., New York, 
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For the thin, elastic, rectangular laterally loaded plate that supported 
the corners, with without tributary beams disposed the panel boundaries, 
various important cases may considered corresponding the limit values 
the flexural and torsional rigidities the beams. Thus, the clamped plate 
supported rigid beams may considered corresponding values the 
flexural and torsional rigidities for the beams that approach infinity. Similarly, 
other combinations limit values these rigidities, which correspond the 
solutions that are known for rectangular plates, are classified Table which 


TABLE DEFINED THE FLEXURAL AND 
RIGIDITIES THE SUPPORTING BEAMS 


Flexural 
(Plain support) (Flexural support) support) 
nonflexural, simple 


nonflexural, torsional 


plain, clamped flexural, clamped nonflexural, clamped 


(Clamped support) 


applies polygons general. Therefore, solution for the laterally loaded 
plate, involving formulas expressing dependence parameters corresponding 
the torsional and flexural rigidities the beams, may considered general 
solution, since would reduce the known solutions for limit values the 
parameters. The solution would also yield classes solutions for intervening 
ranges parametric values. The derivation and application such solution, 
with particular reference the rectangular plate, will investigated this 
paper. 

Historical extensive and distinguished literature the theory 
plates can acknowledged only inadequately the brief list references 
that are limited few distinct solutions, exclusive difference equation 
methods, which have been applied majority the historic problems 
thin, rectangular, transversely loaded plates. 

The first these solutions, for the deflection panels simple, nonflexural 
double infinite series trigonometric terms. Other investigators have applied 
the same form solution the clamped plate, with without nonflexural 
supports. 

The next solution appear, which was also important the sense 
initiating new class applications, was that 1899. The 
latter’s primitive for the applicable differential equation may expressed 


des Recherches sur Flexion des Plans Navier, Bulletin des 
Sciences, Société Philomathique Paris, Vol. 10, 1823, pp. 94-102. 

ibre élastique Plaque Lévy, Comptes Rendus, Vol. 129, 1889, 
pp. 


(Simple support) 
finite plain, torsional flexural, torsional 
(Torsional support) 
m 


ELASTIC SUPPORTS 


the summation; Y,,, the complementary function, comprises hyperbolic func- 
tions the variable the coordinate and particular integral. 
Two applications Lévy’s result, and are also ex- 
amples construction which often recurs the literature. These applica- 
tions involve the superposition pair expressions, each the same form, 
but with and and and mutually interchanged, which and are the 
panel span lengths, measured parallel the and axes, respectively. Mr. 
Hencky employed simple superposition functions (z, and (y, 
which comprised Lévy’s primitive, for the problem the clamped, 
nonflexurally supported panel. Mr. Nadai’s method consisted the construc- 
tion orthonormalized set functions out linear combinations 
comprise the form Lévy’s solution. 

different mathematical method (that the variational calculus) 
third method was introduced Ritz.* this method, series functions 
was set which satisfied the boundary conditions, and the coefficients these 
functions were derived minimum energy considerations. When, however, 
this investigation, the difficulty lies the formulation expressions that 
satisfy the boundary conditions, this method does not appear advantageous. 

satisfy the range boundary conditions considered this investigation, 
sufficient use any convenient primitive the differential equation which 
contains arbitrary constant coefficients. the known solutions, the primitive 
Lévy available for this purpose and may conveniently utilized. 


FORMULATION THE SOLUTION 


the formulation the solution the Kirchoff-Kelvin-Tait theory bound- 
ary conditions has been that, compared with more 
exact assumptions, the possible error this theory small for thin plates and 
the error confined the narrow layers adjacent the boundaries. this 
theory two equations for each boundary edge are required order define the 
boundary conditions for polygonal plate, equations for polygon 
sides. Furthermore, the plate assumed thin, have small deflections, 
and the middle surface assumed inextensional, that sides elements 
that are normal before bending remain normal after bending. agreement 
with these assumptions, the differential equation used the fourth order bi- 
harmonic plate equation 


which the flexural rigidity the plate and the the dis- 


“Der Spannungszustand rechteckigen Hencky, Oldenbourg, Munich, Germany, 

Elastischen Platten,” Nadai, Springer, Berlin, Germany, 1925, pp. 180-184. 

eine neue Methode zur Losung gewisser Variationsprobleme der mathematischen Physik,” 
Ritz, Crelles Journal, Vol. 135, 1909, pp. 1-61. 

Natural Philosophy,” Lord Kelvin (W. Thomson) and Tait, University Press, 
Cambridge, England, Part 11, 1883, pp. 188-192. 

Contributions Applied Edwards, Ann Arbor, Mich., 1949. 
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tributed load. Since the complementary function solution the homo- 


geneous portion this equation includes four constants integration, 


independent solutions superposition are required yield sufficient constants 
for application the equations, namely the case 

Therefore both and are solutions Eq. but may considered 
resulting from rotation Cartesian axes from those w;. Since Eq. 
invariant differential equation under rotation axes, the same form 
when referred the new axes. More explicitly, denoting the new coordi- 
nates primes, and with rotation through 


—the negative signs being included the arbitrary constants implicitly con- 
tained the last equation. 

For convenience, expressions similar and which occur the same 
form, but with both, mutually interchanged, may termed 
being analogous each other. 


The process superposition component solutions referred successively 
rotated Cartesian axes, may extended the construction formal solu- 
tion for the general polygor. However, the contained this paper 
confined symmetrical conditions the rectangular plate. The boundaries 
and orientation the plate are shown Fig. and the positive sense the 
stress resultants shown Fig. flexural rigidity the plate may 

the plate material; the thickness the plate; and denotes Poisson’s 
ratio for the plate material. 
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complementary function Eq. utilized component solution 
Eq. the primitive Lévy: 


42 


which the index being assumed include all integers. Ac- 
cordingly, application Eq. and its analogous expression the expression 
for yields, for symmetrical conditions the rectangular plate under any 
symmetrical loading, 


which the terms containing the arbitrary constants and 
are unaffected the differential operator, and are added for greater 
and Aim, Bim, and are arbitrary constant coefficients corresponding 
each value the index the series comprising the primitive the plate 
differential equation. 

Eq. for each value the index contains four arbitrary constants for 
application the equations which define boundary conditions. will 
shown subsequently, all the classifications Table may expressed four 
equations, for either the symmetrical antisymmetrical conditions the 
rectangular plate. Accordingly, Eq. the solution embracing all the classi- 
fications the boundary conditions listed Table subject the determina- 
tion the corresponding four sets constants, Aim, Bom, and Bim, con- 
sistent with the boundary conditions assigned. 

evident that equations similar Eq. can set for the antisym- 
metrical and also for the general nonsymmetrical plate conditions. 

The set four equations, which represents all the classifications Table 
expresses the equilibrium the plate its boundary edges. These equations 
may written 


Eqs. and are statements for the equilibrium the transverse reaction 
between the plate and the beams the edges and +a/2, 
respectively. The parameters and denote the ratios the flexural rigidi- 


ties the beams the plate the edges and +a/2; that 


Plates and Timoshenko, McGraw-Hill Book Co., New York., Y., 
DP. SY. 
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ticity, and and the moments inertia, the corresponding beams. 
assumed that only forces transverse the plane the plate are transmitted 
between the plate and beams. 

Eqs. and are expressions for the equilibrium the moment between the 
plate and the beams about the axes coinciding with the edges 
+a/2, respectively. The parameters and denote the ratios the 
torsional rigidities the beams the flexural rigidity the plate the edges 
that longitudinal stresses caused warping restraint may neglected the 
beams under torsion. 

The various classifications Table correspond setting zero allow- 
ing the values and and both pairs parameters approach 
infinity. The central position the table corresponds all the finite values 
for all these parameters. 

are the means for the determination the constant coefficients 
for the solution the problem. The process could used straightforward 
manner means obtaining solution having general application all the 
classifications Table However, simplify the presentation, the method 
applied the simply, flexurally supported panel, corresponding line 
Table 


Loap 


The formal solution for the simply, flexurally supported rectangular plates 
under uniform load the inclusive solution given Eq. The particular 
integral utilized, corresponding the uniform loading assigned, 


which s(m) convenient method representing the alternating sign; that 
reduces one determining the values the arbitrary constants contained 
ingly, Eq. becomes for simple suport, 


Substitution 10, which equation for the derivatives Eq. yields 
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Eq. can regarded sum infinite set equations, one equation 
for each value and each equation containing, all its terms, factor com- 
prising trigonometric function similar index. elimination this com- 
mon factor, each equation solved for the corresponding constant, which 
results 


s(m) 


The evaluation expressed Eq. 12a and the trivial values 
are sufficient satisfy Eq. for all values Thus the monomial 
terms containing and are eliminated from the solution for problems 
simple support. 
Similarly, the equation analogous Eq. 12a, obtained application 


1 — 


Eq. becomes, substitution the appropriate derivatives 


m~1,3,°-- b 


aa 


m=0,1,2, 
m=0,1,2,°°° 


The series terms Eq. 13, comprising mixed hyperbolic, trigonometric, 
and other functions, not enable the solving for the unknown constants 
direct manner. However, process utilized among others, the 
mixed terms are first expressed equivalent trigonometric series, with period 
2a, which results 


which s(n) and the indices are assumed com- 
prise all integers, except otherwise indicated. 
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The order summation the double series shown Eqs. may re- 
versed interchanging the indices and substitution Eq. 
these series, with interchanged indices, and elimination the factor cos 
the equation for each odd value 


n=1,3, 


For each even value the index the corresponding equation homo- 
geneous, that the right-hand side Eq. vanishes. Hence, sufficient 
satisfy the equations for the boundary conditions the trivial solutions for 
the constants even index; that is, 


Accordingly, the Fourier coefficients and Eqs. 14, with inter- 


changed indices, are evaluated only for odd values and corresponding 
Eq. 15, use the expressions 


Qinm 


m? a? a’ 

These expressions for the Fourier coefficients and Eqs. for Aim and 
are simultaneously introduced into Eq. 15, and the terms resulting coeffi- 


cients are simplified. further simplification the result use the 
summation 
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and the notation 


= 2 3 


and 


the dimensionless form derived from Eq. 


Bom Bin 


equation analogous Eq. 20a can derived from Eq. and the 
form 


The constants Bim, with equal which, with the related 
constants and were required satisfy the equations for the boundary 
conditions assigned, comprise the formal solution the infinite sets numerical 
equations represented 20. 

practice, the unknowns Eqs. are approximated solving, 
conventional manner, partial sets the numerical equations. illustration 
the procedure for the numerical application Eqs. provided the 
following example, that corresponds arbitrarily assigned value the 
parameter 

Numerical simplify the example, the plate assumed 
square, and the sides are equal The parametric value assigned the 
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reduced the same form, which for the numerical values assumed become 


68m 


partial set numerical equations obtained from Eq. 21a, corresponding 
the first six values the index that For load 
these equations become 
—0.8350 
—0.5015 
—0.1289 
—0.5115 


—0.2533 
—0.1435 


This set numerical equations may represented algebraically the 
relationship which the matrix the numerical coefficients the 
left-hand side Eqs. the column matrix comprising the unknowns B,,, 
with 3,---, and ris the column the right-hand side The 
solution for may obtained iterative operational formula first re- 
ported and the form: 


provided the matrix series indicated converges, which the matrix minus 
the unit diagonal; that ist Accordingly, substitution the corre- 
sponding matrices into Eq. 22a yields successively 


These values the unknowns, substitution into Eqs. 21b agree within 
one unit the last significant digit the respective right-hand side. The cor- 


Eqs. and and converge rapidly, follows 


—0.8001 


These numerical values, multiplied the applied load are inserted for 
the corresponding constants into Eq. The result the evaluation 
partial sum for the assigned values the parameter 


Simple Method for Solving Simultaneous Linear Equations Successive Approximation 
Morris, Journal, Royal Aeronautical Soc., Vol. 39, 1935, 349. 
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Corroboration Known can verified that for and 
Eqs. reduce expressions that are known literature for this 


limit case nonflexural support. 
Thus for from Eq. 20a, 


cosh 


and from Eq. 12a 


Similarly, expressions analogous Eqs. and for the constants and 
Aim are derived use Eqs. 20band 12b. These expressions for the values 
Bom, Bim, Aom, Aim, together with the formal expression for Eq. 
complete the required solution. However, for the nonflexural case, Eq. may 
simplified since each solution for and separately satisfies the boundary 
conditions all four sides. Since there can only unique solution, 


m=1,3,-°° 


which the constants and are evaluated Eqs. 23b and 23c and Aim 
and Bim the analogous expressions. Thus, the left-hand side Eq. 24, 
multiplied two, agrees with classical solution for this case nonflexural 
support." 

can shown that for the classification and finite, Eqs. 
also reduce corresponding expressions derived the literature. 

Other Assigned other load functions, particular integrals 
expressed infinite trigonometric series may introduced for the formal 
solution, and equations similar Eqs. can obtained the same method 
derivation. The left-hand sides Eqs. remain the same for all particular 
integrals expressed infinite trigonometric series, evaluated for odd integer 
values the index. 


Plates and Timoshenko, Book Co., Ine., New York, Y., 
Pp. 
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NUMERICAL EVALUATIONS AND REACTIONS 


Numerical evaluations were made the author various moment and 
reaction coefficients for square plates over the full range relative rigidities 
the supporting beams. Rectangular plates, simply supported and uniformly 
loaded, were also investigated. 

For square plates, values the parameter were selected intervals 
the range from and corresponding sets computations the 
constants B,, and were carried through described previously. However, 
for greater accuracy, summations were carried for ten values the index, and 
sets ten numerical equations each case were solved for the constants B,, 
with 19, method successive elimination adapted for 
computing machines Crout.” 

Sets moment and reaction coefficients were computed for ten distinct 
stations along the axes symmetry, along the diagonals, and along the edges 
the panel. The differential expressions used for the computation the mo- 
ments comprise the left-hand side Similarly, expressed 
the left-hand side The differential expressions used for the reaction 
the edges b/2 and a/2 comprise, respectively, the left-hand sides 
Eqs. and The evaluation for each station the derivatives ap- 
pearing the differential expressions, was effected substitution the con- 
stants A,, and B,, the respective expressions for the derivatives and 

The results these computations are shown Table and Figs. 


0.25 0.36125 


0.0635 0.0317 0.0187 —0.02834 —0.03941 
0.026340 0.020587 0.015981 0.013021 0.009044 0.006015 0.004062 


P at/N 


ry 


For square, simply supported plates under uniform load Moment coefficients 


refer twisting moments the corners the panel and the deflections are the center the panel. 


Particular panels corresponding the parametric values 


found that the various expressions for moment reaction reduce simple 
exact formulas. Thus, all cases which correspond Eqs. 25, the right-hand 


Short Method for Evaluating Determinants and Solving Systems Linear Equations with Real 
Crout, Transactions, Am. Inst. Electrical Engrs., Vol. 60, 1941, pp. 
1235-1240. 


0.0 0.10 1.0 3.0 
and 


0.25 
Values 

0.25 
Values 


4.—Moment 


0.20 


0.15 


Positions 


Coordinate 
Corresponding Coordinate 
Positions 


god 


0.16 
0.14 
0.12 


Corresponding Coordinate 


Positions 


od 


Values 


2 
a 


od 


(DIAGONAL) 


0.02 
Values 
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Reaction Coefficients C,(Edge) 
Values 


Values 


sides Eqs. vanish, and hence all these cases the trivial values 

Bim with are sufficient satisfy the equations for 

the boundary conditions. For all these cases, the corresponding constants 

obtained are 


sech 


The terms containing which appear the various expressions 
for moment and reaction are further simplified for the values indicated 
26, with the aid the relationship 


0.36125 
020 025 030 035 040 045 0.50 
and 
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The derivation Eq. shown Appendix together with the applica- 
tion this equation reducing the expressions for moment and reaction the 
following simple formulas (corresponding Eqs. 25). 


m~1,3, 
and 


2 
Analogously, 


Eqs. can corroborated the procedure previously described 
for sets numerical evaluations, with numerical values determined from Eqs. 
substituted for the constants polynomial expression for conform- 
ing Eqs. 25, and Eqs. 31, constructed independently the general 
solution, provided Appendix 


EXPERIMENTAL VERIFICATION 


Experimental verification the analysis previously presented was sought 
two introductory tests plates simply supported beams arbitrary 
flexural rigidity. Glass and gypsum plaster were selected plate materials, 
because the ideal linearity the stress-strain relationship the point 
rupture for both materials. was thought that the test results could sim- 
plified the use plaster that would allow for the casting the model with 
the supporting beams embedded symmetrically within the slab thickness. 
Also, from preliminary tests with sample specimens, was concluded that com- 
merical resistance type strain gages are convenient use with plaster models and 
consistent results. 

Model Specifications and Test test model dimensions and 
the details support are shown Fig. The uniform loading was applied 
each test water contained rubberized fabric tank, placed contact with 
the plate. The tank portions extending beyond the panel boundaries were sup- 
ported independently the panel order eliminate any vertical component 
pull the panel edges. Strains the plate and beam surfaces were ob- 
served means SR-4 resistance strain gages, type A-5, bonded the plates 
the positions shown the plan views. 


2 


ELASTIC SUPPORTS 


Test plate material consisted commercial plate glass, nominally 
}-in. thick. The four supporting bars were steel, with mitre cut the ends 
and arranged frame, shown Fig. 7(b). The supports for the bars 
each corner the frame consisted small disc, which was laid flat rigid 
plane base. For this preliminary test, the gages were applied one side the 
glass plate and the supporting bar. The glass was first placed with gages 


(a) PLAN VIEW, TEST PANEL (c) PLAN VIEW, TEST PANEL 


Flanges 
Truss Shaped Open Web 


Plate Glass Gypsum Plaster 
Average Slab 0.472 In. Average 


Thickness Thickness 


Center Lines Beam 


Supports Beam Ends 


SECTION, TEST PANEL (d) SECTION, TEST PANEL 


the compression side. After taking initial readings, strain observations were 
recorded corresponding three increments water head. The water was 
then removed, the plate reversed, successive increments load reapplied, and 
corresponding tensile strains recorded. 

Test gypsum plaster used for the plate material was commercial 
variety with water-plaster ratio 70:100 weight. The supporting 
beams, which were fabricated with open webs, were embedded the cast slab 
along the edges, and were connected the panel corners small blocks resting 
jacks. Strain gages were applied the top and bottom surfaces the panel 
and supporting beams, usually pairs each gage station. Other gages, 
not indicated the figure, were applied the panel close the beams, and 
the beams close the ends, determine whether simple support was 
realized. 

Successive increments water head were applied and corresponding strain 
readings recorded. Following each increment head, the water was entirely 
removed ascertain whether the strain readings reverted the initial zero 
values with each release the load. 


23.4 In. 23.0 In. 
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The modulus elasticity the plaster was determined specimens the 
shape small flat beams cast from the same batch used for the test panel, 
and also specimen cut out from the test panel the conclusion the ex- 
periment. From tests these specimens flexure, mean value was deter- 
mined for the modulus elasticity the plaster 1.22 10° per in. 
Poisson’s ratio was determined experimentally 0.14. 

Test Model Values the glass Test the modulus elasticity 
was assumed per in. and Poisson’s ratio was assumed 
0.21. The modulus elasticity for steel 10° per in. Accord- 
ingly, for glass plate 0.240 in. thick, 0.00120 For the supporting bar, 
0.00200 and A/a 0.214. 

similar way, with the aid the experimentally derived values and 
for plaster, the relative flexural rigidity the plaster model was \/a 0.407. 

Evaluation Comparative Test strain readings, observed 
the laboratory, were reduced equivalent mean strain values per unit 
applied load, for uniform thickness plate each test. The observed strain 
values were corrected for variation actual thickness multiplication 


the original readings the factor the case glass plate, and 


0.472 

for the plaster plate. These corrected strain values were reduced 
strain per unit pound per square inch application numerical factors cor- 
responding each increment water head. Finally, for each gage station, 
the algebraic mean was obtained for the corresponding compression and 
tensile (+) strains, corrected for thickness and for unit load. These mean 
strain values comprise the final experimentally derived strains each gage 
station, which were plotted, for comparison with the theoretically derived 
strains, Fig. 

The algebraic mean values strain thus obtained part eliminate 
T-beam effect, whereby the strains the top surface the plate were con- 
sistently higher than the strains the bottom surface, the same gage 
stations. This effect was found during both tests and attributed the 
design the panel supporting members. Accordingly, the test results shown 
Fig. include only algebraic mean values derived from gage readings where 
available reading for the top and one reading for the bottom sur- 
face, respectively, plate beam each gage station. 

For Test the values plotted are the averages these algebaric mean 
values for corresponding positions where available two orthogonal axes, for 
various intensities load. These average strain values further eliminate 
effects which are unsymmetrical with the two axes. 


Theoretically Derived Strain value 0.214, evaluated 


means the dimensions and elastic constants model agrees closely with 


for which value the parameter numerical example was solved. 


Hence, corresponding 0.2125, theoretical strains normal 
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axis the plate, were computed for points generally along the axis. These 
computations can made with the aid moment coefficients computed for 
corresponding positions, means the relationship 


a? 
2 oy = ON (Cow — # oe (32) 


which substitution made for the dimensions and the ratio conform 


the test model. Theoretically derived strains thus computed for the condi- 
tions Test are plotted Fig. for comparison with test results. 


Strain Microinches per In. 
per per In. Applied Load 


Test Experimental Strains One Axis 
Average Test Strains Two Axis 
Test Strains One Axis 


Values 


Fra. 9.—A Comparison Test AND THEORETICAL STRAINS 


Similarly, 0.407, for model approximately equal the parameter 
value 0.36125, for which the moment coefficients are expressed exact formu- 
las given Eqs. and 29. Hence, the moment coefficients for 0.407 


were obtained use Eqs. and 29, and corrected interpolation the 
sets coefficients plotted Figs. The strains for the conditions 
Test were then computed means Eq. and plotted Fig. 

Comparison Between Test and Theoretical Values.—For Test the average 
deviation the test strains plotted Fig. with respect the corresponding 
theoretical values represented the curve, evaluated 12.6%, based 
the theoretical results. 

For Test the average deviation between theoretical and test 
values 9.1%, where these test values represent averages derived from gage 
readings corresponding positions two orthogonal axes the model. How- 
ever, the average deviation thus computed does not include considerable devi- 
ation, amounting 23.6% for position the test model where only uniaxial 
strains were available. this position, the large deviation represents, part, 
apparent error due the unsymmetrical strains experienced model 


500 
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The curves Fig. refer only points along the axes parallel the panel 
edges. For the station the diagonal axis, shown the panel diagram 


0.8 (approximately the quarter point) the theoretical and experimental 


values are, respectively, 476 and 466 microinches per in. per per in. 
applied load, where however, the experimental values are based only uni- 
axial gage readings. 

Test Conclusions.—The average deviations evaluated correspond average 
test strain values which are 87% 91% the respective theoretical values. 
This numerical comparison, together with the graphical comparison Fig. 
leads the conclusion substantial corroboration the analysis. 

Moreover, the deviations found between theoretical and experimental strains 
are attributable, part, the discrepancy between the simplified boundary 
conditions assumed the analysis (of vertical forces between plate and beams) 
and the T-beam condition actually experienced. concluded that the 
theoretical strain values computed according the methods presented herein 
approximate the algebraic mean values compressive and tensile strains de- 
veloped the corresponding surfaces panels with T-beam action. 

future tests, for the sake greater accuracy correlation between test 
and theoretical results, further precautions shou!d taken for elimination 
T-beam action. 


CONCLUSIONS 


This paper presents method for the solution boundary value problems, 
which adapted problems transversely loaded plates, well other 
physical problems, especially those involving harmonic differ- 
ential equations. 

The method applied the formulation solution for the problem 
the rectangular supported beams along all the edges, for general values 
parameters expressing the flexural and torsional rigidities these beams. 
Most the solutions the this subject apply limit values 
these parameters. far known, the solution this paper has not been 
obtained other investigators. 

The formal solution comprises infinite series terms with constant co- 
efficients, dependent the boundary conditions, most solutions plate 
problems literature. However, the determination these constants involves 
the solution sets simultaneous linear equations. These constants, once 
obtained for the assigned values the parameters, and for the assigned plate 
serve thereafter for computation stress and deflection coeffici- 
ents for various positions the corresponding panels. found that the 
convergence the constants and with 1,3,--- very rapid. 
These constants, used expressions for shear and moment, occur com- 
binations terms, the convergence which, extreme cases, isslow. these 
extreme cases, however, the accuracy still estimated adequate for design 
graphical purposes. Similar conclusions accuracy apply the limit 
case, which the solution for which identical with the classic solution 
for nonflexural supports. 
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The detailed application the formal solution, complete with numerical 
results corresponds the problem simple, flexural support under uniform 
load. For this problem, comparative numerical investigation was conducted 
moment and reaction coefficients relation flexural rigidity supporting 
beams. The graphs clearly show the considerable moments developed along 
the edges which comprise characteristic flexural support. Thus, 

mum along the edges. The pattern reaction along the edges 
changes from downward convex downward the parameter lambda 

equal its value give rise progressively larger changes moment 
and shear. 


all square panels having relative the moments are maxi- 


For any rectangular panel, corresponding the particular values 


exact formulas that represent parabolic cylinders superposed symmetrically 
with the axes symmetry the 

Maximum defllections, which are general interest, are among the quanti- 
ties that are easily obtained with extreme precision. From Table these 


quantities are seen vary from 0.026340 for 0.004062 for 

Introductory experimental results are found consistent with the the- 
oretically computed corresponding data. The test results obtained also sug- 
gest improvements the preparation and arrangement panel specimens 
future tests, provide for increased accuracy experimental corrobora- 
tion the theory. 


the derived expressions for moment reduce simple 
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APPENDIX DERIVATION EQUATIONS 


Since application the same differential operator and 
yields 


Substitution Eq. 33, termwise differentiation Eq. 24, results Eq. 27. 
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The following equations are obtained from Eq. termwise differentia- 
tion integration; that is, the differentiation Eq. with respect 


(m) 


m~1,3,++- 


and the integration Eq. with respect 


m~1,3, 


and the differentiation Eq. 34b with respect 


m~1,3, 


may noted that for the constants Bim and are the equiva- 
lent, except for constant factor, Eqs. for the constants and 
Hence, substitution into the differential expressions for moment, 


M, N (Wyy + Wz2) 


found, with the aid Eq. 27, that all the terms containing and 
vanish from the resulting expressions, which readily reduce Eqs. and 29. 

Similarly, Eq. 31, for twisting moment, readily obtained with the aid 
Kq. 34c, from the differential expression M,, 

Also, Eqs. are readily obtained with the aid Eq. 34a from the differen- 
tial expressions represented the left-hand sides Eqs. and 

The polynomial expression for constructed satisfy the boundary con- 
dition (Eq. 10) for simple support: 
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found that Eqs. and for equilibrium are also satisfied Eq. 35c, 
provided that the values for and conform Eqs. 25. fol- 
low from Eq. 


APPENDIX II. NOTATION 


The following letter symbols, adopted for use the paper and for the guid- 
ance discussers, conform essentially with Standard Letter Sym- 
bols for Structural (ASA Z10.8—1949), prepared committee 
American Standards Association with ASCE participation, and approved 
the Association 1949. 


Aim, Aom, Bim, Bom arbitrary constant coefficients corresponding each 
value the index the series comprising the primitive the plate 
differential equation. 

panel span length measured parallel the z-axis. 
panel span length measured parallel the y-axis. 
moment coefficient for the moment M,. 
moment coefficient for the moment 
reaction coefficient. 
column matrix comprising the unknowns B,, with 
modulus elasticity the plate material. 


modulus elasticity supporting beam material edges 


depth thickness the plate material. 
unit diagonal matrix. 
moment inertia the cross-sectional areas the supporting 
moment inertia the cross-sectional areas the supporting 


ko, constants used the general equation. 
matrix the numerical coefficients the left-hand side Eq. 21b. 

bending moment per unit length about axis parallel the y-axis, con- 
sidered positive when causing compression the top fibers 
the plate beams. 

bending moment per unit length about axis parallel the z-axis, 
same sign convention for M,. 

twisting moment per unit length section plate normal the 
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twisting moment per unit length section plate normal the 
y-axis. 
index summation. 
flexural rigidity the plate material. 
index summation. 
intensity the distributed load, positive downward. 
right-hand side Eq. 
matrix minus the unit diagonal 
deflection the plate, positive downward. 


strain the plate. 


ratio the torsional rigidity the flexural rigidity edge 


a 


Poission’s ratio. 
particular integral. 
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